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O ' ABSTRACT 

(N ! 

High redshift radio galaxies are great cosmological tools for pinpointing the most massive objects in the early 
Universe: massive forming galaxies, active super-massive black holes and proto-clusters. We report on deep 
narrow-band imaging and spectroscopic observations of several z > 2 radio galaxy fields to investigate the 
nature of giant Ly-a nebulae centered on the galaxies and to search for over-dense regions around them. We 
discuss the possible implications for our understanding of the formation and evolution of massive galaxies and 
> ■ galaxy clusters. 
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1. RADIO SOURCES AS COSMOLOGICAL PROBES 

High redshift radio galaxies (HzRGs; z > 2) are great beacons for pinpointing the most massive objects in the 
early universe, whether these are galaxies, black holes or even clusters of galaxies. At low redshifts powerful, non- 
thermal radio sources are uniqely associated with massive ellipticals. Their twin-jet, double-lobe morphologies 
and large luminosities suggested already early on that such galaxies must also have spinning, super-massive 
black holes (SMBH's) in their centers. 1 ' 2 We now know that the masses of the stellar bulges of galaxies and 
their central black holes are correlated, 3-5 suggesting a causal connection. If radio sources are powered by 
SMBH's then it is no longer a surprise that their parent galaxies occupy the upper end of the galaxy mass 
function. 

There is excellent evidence that radio galaxies are also the most massive systems at high redshifts, even 
though their parent galaxies are very young and may still be forming. The combined near-infrared 'Hubble' 
K — z relation for radio and field galaxies 6 shows that HzRGs are among the most luminous systems at any 
given epoch up to z ~ 5. Between < z < 2.5 this K — z diagram can be modeled using passive evolution of 
5£* ellipticals 7 which formed at 5 < z < 10. Other evidence includes the direct detection of absorption lines 
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from massive young stars, 8 large (~ 50 — 70 kpc), multi-component rest-frame UV and optical morphologies, 9 
the 'hyper' luminous rest-frame far-infrared luminosities (Lfir 10 13 Lq) and huge implied star formation 
rates 10 ' 11 (> 2000 M Q yr -1 ), and last but not least, the very extended (~ 30 — 50 kpc) molecular gas and dust 
clouds that have recently been discovered 12 ' 13 around several HzRGs showing that the star formation occurs 
on galaxy wide scales. 

HzRGs are also excellent tools for finding over-dense regions ('proto-clusters') at high redshift. This is 
because, in standard Cold Dark Matter (CDM) scenarios, galaxy formation is a highly 'biased' process: the 
most massive galaxies, and the largest clusters of galaxies, are expected to emerge from regions with the largest 
over-densities. 14 Simply put: the most massive systems (galaxies, SMBH's and galaxy clusters) hang out 
together, and radio sources are a great way to find them and to investigate their interrelations and evolution. 

In this paper we report on deep narrow-band imaging and spectroscopic observations of several z > 2 radio 
galaxy fields using the Keck and ESO/VLT telescopes to investigate the nature of giant Ly-a nebulae centered 
on the galaxies and to search for over-dense regions around them. We discuss the possible implications for 
our understanding of the formation and evolution of massive galaxies and galaxy clusters. We will adopt the 
cosmological parameters Om = 0.3, Qa = 0.7, H = 65 km s -1 Mpc -1 , for which the age of the Universe at z ~ 
2, 3 and 4 is 3.5, 2.3 and 1.6 Gyr respectively, and the angular-to-linear transformations are 9.0, 8.3 and 7.5 
kpc arcsec -1 . 

2. GIANT Lya NEBULAE 

Since the discovery of the first HzRGs it has been known that they are surrounded by large Lya emission line 
nebulae. 15-18 Because the nebulae are centered on large, young galaxies they provide a unique opportunity 
to study how such galaxies form (accretion or merging ?), and about AGN/starburst feedback and chemical 
enrichment during this process (heating or outflow ?). We therefore observed the Lya nebulae of several HzRGs 
and describe the results for two of them in detail below. A summary of the largest high redshift Lya nebulae 
currently known is given in Table 1. 

Table 1. Large, High Redshift Lya Nebulae 
Name z Size Log(Ly-a) 6 Telescope Reference 
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MRC 1138-262 


2.156 


150 


X 


120 


44.3 


VLT 8-m 


Rcf. 


18 


Francis' Blob" 


2.380 


100 


X 


30 


43.9 


CTIO 4-m 


Ref. 


19 


Steidel Blobs" 


3.09 


140 


X 


120 


44.1 


Palomar 5-m 


Rcf. 


20 


B2 0902+34 


3.395 


140 


X 


100 


44.8 


Keck 10-m 


Rcf. 


11 


4C 1243+036° 


3.570 


160 


X 


50 


44.0 


ESO 4-m 


Rcf. 


17 


4C60.07 


3.791 


110 


X 


60 


45.1 


Keck 10-m 


Rcf. 


11 


4C41.17 


3.798 


190 


X 


130 


45.2 


Keck 10-m 


Ref. 


11 


TN J1338-1942 


4.102 


150 


X 


40 


44.7 


VLT 8-m 


Ref. 


21 



a Radio 'quiet' 

6 Depends on surface brightness limit (larger telescopes detect larger 
Ly-a nebulae) 

c Maximum size determined from long slit spectrocopy, not imaging 



2.1. 4C41.17 at z = 3.798 

4C41.17 was imaged using a custom-made, high throughput interference filter with a 65 A bandpass centered 
at the redshifted Lya = 5839 A line with the Echcllcttc Spectrograph and Imager 22 at the Cassegrain focus of 
the Keck II 10m telescope. The data were obtained during photometric conditions and good seeing (FWHM = 
0.57"). With a total exposure time of 27,600s this Lya image is the most sensitive obtained to date and reaches 
a surface brightness of 8.0 x 10 -19 ergs -1 cm -2 arcsec -2 (3<r limit in a 2.0" diameter aperture). A broad-band 
R image, obtained as part of a multi-band spectral energy distribution study of the 4C41.17 field was scaled 




Figure 1. [Left] Keck narrow-band Lya image (grey scale) of 4C41.17 overlaid with a 6 cm VLA radio image (contours). 11 
The cross shows the position of the radio core. [Right] Normalized surface brightness profiles of the [Oil] (solid) and Lya 
(dashed) emission along the inner radio axis and the long south-west filament. Top: Extended [Oil] in the velocity range 
to +500 kms -1 . The spatial zero-point corresponds to the position of the radio core. Middle: Similar to top panel but 
for the velocity range —500 to — 1200 kms -1 (Note: The to — 500 kms -1 range for [Oil] is affected by near-infrared 
sky lines and is not shown). Bottom: Relative velocity (solid line) and velocity dispersion (bars) of the Lya emission. 
Bar and symbol size indicate the respective uncertainties in the individual measurements. The projected distances of 
the south-west and north-east radio lobes along the slit direction are indicated (dotted lines) . 



and subtracted from the Lya image to construct the pure emission-line image shown in Figure 1. With linear 
dimensions of 190 kpc x 130 kpc at the detection limit the 4C41.17 nebula is the largest presently known (Table 
!)• 

A striking morphological feature of the nebula is a cone-shaped structure emanating from the center of the 
galaxy with a 75 kpc long radial filament in the vicinity of the south-west radio lobe and a crescent-shaped 
cloud with radial horns. This morphology resembles that of other, nearby active and starburst galaxies and is 
suggestive of entrainment by the radio source, 23 outflow driven by radiative pressure from the AGN, 24 or a 
starburst superwind. 25 

Since Lya is a resonance line, it is important to determine whether the nebula is ionized, or rather, results 
from the scattering of the Lya photons (produced near the nucleus) by an extended neutral hydrogen gaseous 
halo. We therefore obtained a long-slit spectrum of 4C41.17 with the Keck II near-infrared spectrograph 26 
(NIRSPEC) to measure the extent and kinematics of other lines, [O II] A 3727 and [O III] A 5007, along the 
brighcst filament of the Lya nebula. We discovered both extended [Oil] and [OIII] emission. In particular, 
[Oil] emission was detected as far as <~ 60 kpc from the nucleus (Fig. 1). This shows that the Lya nebula, at 
least in this direction, is not due to scattering but must be locally ionized. Furthermore, the presence of the 
emission lines of oxygen shows that the halo gas is not chemically pristine (primordial H, He) gas but has been 
enriched either at much earlier epochs (during the 'Dark Ages'), or by more recent star formation. 

The near-infrared [Oil] and optical Lya kinematics show that both emission lines exhibit large blue-shifted 
velocities ~ 600 — 900 km s -1 (in projection) along the radio axis. In particular, the gas is very disturbed along 
the south-west filament, with Lya velocity widths ranging up to Avf ~ 900 — 1600 km s -1 . Beyond the radio 
hotspot, the velocity and velocity widths decrease abruptly. Evidence for entrainment of emission-line gas in 
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Figure 2. [Left] Keck narrow-band Lya image 11 (grey scale) of 4C60.07 overlaid with a VLA radio image (contours) at 
6 cm (contours). The cross shows the position of the radio core. [Right] HST 'R'-band image of 4C60.07 (grey scale) 
overlaid with a Keck narrow-band Lya image (contours) and IRAM interfermometer dust image 11 (dotted lines). 



long filaments, even if at considerable distance from the observed radio emission, has also been seen in nearby 
radio galaxies i.e., 4C29.30. 23 

We note that in the canonical picture of radio sources the hotspots and lobes are surrounded by bowshocks 
and cocoons of radio quiet, shock heated gas with a scale sizes wich are significantly larger than the observed 
radio emission. 27, 28 In this picture the emission-line filaments located at the interface of the cocoons and 
the ambient gas and are not in direct contact with the radio lobes or hotspots and can even extend beyond 
radio hotspots 29 (Figs 1, 2), in particular if one also accounts for projection effects and the fact that the radio 
observations only show the highest surface brightness regions. The kinematics along the filament in 4C41.17, 
its radial filamentary structure, and the chemical enrichment of this gas therefore all indicate a process of 
entrainment of material away from the central regions by the radio source. The outer regions of Lya nebulae 
may not have been affected by the sources, especially if they are still young and small, and one would expect 
these regions to be less enriched, as has been observed. 30,31 

2.2. 4C60.07 at z = 3.790 

Because of its very similar redshift the radio galaxy 4C60.07 could be observed with the same narrow-band filter 
as 4C41.17. The total observing time on source was 7,200s, resulting in a surface brightness detection limit 
of ~ 2 x 10~ 18 erg s" 1 cm -2 arcsec -2 (3<r limit in a 2.0" diameter aperture). Figure 2 shows the Lya nebula 
in relationship to the radio source, rest-frame UV HST structure, and the extended dust and molecular gas 
that has been found in this system. 11, 12 Unfortunately the Lya image is not as deep as for 4C41.17, not only 
because of the shorter intergration time, but also because 4C60.07 has significant Galactic foreground extinction 
(~ 1.6 magnitude at the central wavelength of the narrow band filter). 

Nevertheless one notes several similarities to 4C41.17 and other radio galaxy emission-line nebulae. First, 
the brightest Lya emission is associated with the radio lobe which is closest to the AGN. This is presumably 
because of asymmetries in gas density distribution, with denser gas providing more rampressure on the hotspots 
of the expanding radio source. 32 Second, the galaxy structure is very clumpy and approximately aligned with 
the radio source. This may be due to induced star formation by radio lobes which expand sideways into the 
cold, dense gas seen in the mm-maps. Good evidence for this has been found in 4C41.17. 8,33 Third, and 



this is the first such evidence, the Lya and dust/molecular gas emission show little overlap. We don't know 
how common this is since to date only three HzRGs have been imaged with mm interferometers with sufficient 
sensitivity to map their cold gas and dust emission and of these only 4C60.07 has a deep Lya image. However, 
an anti-correlation between Lya and cold gas certainly is not suprising given that Lya is easily scattered and 
absorbed. In fact the overall impression is that the Lya emission is more or less orthogonal to the dust and 
cold gas distribution, resembling M82 and its emission-line starburst superwind. 25 

2.3. Conclusions for Ly— a nebulae 

Since many HzRGs seem to be surrounded by giant Lya nebulae it seems likely that they play an important 
role in the formation of massive galaxies. What is the origin of this gas, what its fate? 

With respect to the origin we should note that large Lya nebulae also exist without central radio galaxies or 
AGN, and that these also show large velocity widths and gradients. 19,20 In the absence of an obvious central 
source of ionization and/or outflow (but see Ref. 34) this suggests that Lya nebulae are due to the accretion 
of primordial (?), cooling gas in large CDM halos. 35 ' 36 This gas provides a large reservoir of galaxies which 
supplies the building materials for galaxies forming at their centers. 

As to the fate: Starburst 37 and AGN 38 outflows may eject accreting gas in forming galaxies when such 
outflows can overcome the galaxian potential, thereby limiting star formation and galaxy growth. This feedback 
mechanism may explain the fixed black hole to bulge mass ratios observed in nearby galaxies. In the case of 
HzRGs our data shows that this outflow may be significantly helped by shock heating and entrainment from the 
radio sources. We note that radio source induced 'hot bubbles' have also been invoked to reduce gas accretion 
in the centers of nearby cooling flow clusters, which seems required to explain X-ray observations. 39,40 

Radio sources may be a significant source of heating and chemical enrichment during their life time 41 
(~ 10 6 — 10 7 yr). The cross sections of radio lobes are large 27 : tens of kpc 2 and much larger than the hotspots, 
which are only the highest surface brightness features that are visible in radio maps because of the steep radio 
lobe spectra and high restframe frequencies. The multiple component, asymmetric, and twisted radio structure 
of 4C41.17 and many other HzRGs suggests that the central SMBH may experience multiple periods of radio 
source activity and precession, presumably triggered by the interaction with one of the many clumpy components 
which make up the complex galaxy structure. If individual clumps move at radial velocities comparable to the 
velocity dispersion in the proto-clusters (~ 500 kms -1 , Table 2) and the closest are within ~ 10 kpc of the 
SMBH 9, 11 then it takes only ~ 2 x 10 7 yrs to re-supply the AGN. The duty cycles for radio source activity 
could therefore be very short, at least during the period that the galaxies are being assembled. During this 
period a considerable fraction of the nebular gas would be heated and expelled. Not only could this then be 
the 'end of the beginning' for the galaxy, as its collapse gives way to mass ejection, it would also be the 'end 
of the beginning' for the proto-cluster as more and more metals are dispersed throughout the gas, stimulating 
cooling, star formation and galaxy evolution. 

3. PROTO CLUSTERS 



Table 2. Proto-Clusters 



Name 




z 


N a 
N B 


N b 

spec 


a ( km s 1 ) 


FOV c 


Telescope 


Reference 


MRC 1138- 


-262 


2.156 


51 


15 


280, 520 


35 


VLT 8m 


Ref. 42 


53W002 




2.39 


14 


5 


530 


48 


KPNO 4m 


Ref. 43 


MRC 0943- 


-242 


2.919 


55 


18 


824 


40 


VLT + Keck 


Ref. 44 


SSA22a 




3.09 


72 


10 


400 


77 


Palomar 5m 


Ref. 20 


4C41.17 




3.800 


18 


3 




3.5 


Keck 10m 


Ref. 11 


TN J1338- 


1942 


4.102 


28 


20 


325 


40 


VLT 8m 


Ref. 21 



a Natb = Number of Ly-a excess candidates 

6 N spec = Number of spectroscopically confirmed Ly-a galaxies 

c FOV = field of view in sq.arcmin. 



Clusters of galaxies are convenient targets for studies of the formation and evolution of galaxies and may 
help constrain important cosmological parameters. Optical and near-IR surveys, assisted by X-ray selection, 
have found clusters 45, 46 out toz~ 1.3. Very deep optical 'Lyman break' surveys have discovered galaxy over- 
densities ('proto-clusters') at z <~ 3. 09. 20,47 Such observations are time consuming and increasingly difficult 
at higher redshifts. They will not be able to provide substantial numbers of over-dense regions spread over a 
large range in redshift, which would be needed if one wishes to investigate how proto-clusters evolve, at what 
redshifts galaxy clusters become virialized, and what the cluster mass functions are. 

A practical way to identify proto-clusters is to search for over-dense regions around targeted objects. Radio 
sources are ideal for this since they are associated with the most massive galaxies and, in hierarchical CDM 
scenarios, therefore live in the most over-dense regions. Furthermore, radio selection is not biased towards 
dusty systems and provides an important alternative method to the optical/near-IR (rest-frame UV/optical) 
color selected searches. 

The quickest way to identify potential over-dense regions around HzRGs is to use narrow-band filters, or 
better even, 'tunable' filters or other such imaging devices (Section 4) if they were available on the big telescopes. 
At least some fraction of the young galaxies in proto-clusters can be expected to be active starforming systems 
and such Lya excess galaxies (LaEG's) should be detectable in deep emission-line searches (Lya at optical 
wavelengths for z > 2 objects). A deep narrow-band Lya image of the z = 3.09 SSA22a proto-cluster field, 
which was initially found using broad-band color selection, confirmed this. 20 

Our Leiden, LLNL and ANU/AAO groups have undertaken a joint program of deep narrow-band and 
tunable filter imaging of HzRG fields with redshifts ranging between 2.2 < z < 5.2 to identify proto-galaxy 
candidates through excess Lya emission. This is then followed by multi-slit spectroscopy to confirm that their 
redshifts are close to that of the target HzRG. This method has been highly successful and we discuss two recent 
results in more detail below. A summary of presently confirmed proto-clusters is shown in Table 2. 

3.1. MRC 0943-242 at z = 2.919 

MRC 0943—242 was imaged using a narrow-band filter with a 68 A bandpass centered at 4781 A with the 
ESO 8.2m VLT Kueyen telescope using the imaging mode of the FOcal Reducer/low dispersion Spectrograph 
2 (FORS2). The total exposure time was 22,500 s, reaching a magnitude limit of NB^s = 26.7 per (3<r, 
2.0" aperture) and the data were obtained during photometric conditions and 0.8" seeing. A broad-band B 
image (Bab = 27.0, 3a, 2.0" aperture) was used to identify Lya excess galaxies in a manner as described in 
previous papers. 21 ' 42 A total of 77 Lya excess galaxy candidates were found above a limiting flux density of 
5 x 10~ 18 ergs -1 cm~ 2 arcsec~ 2 (3 cr; Fig 3). 

We subsequently obtained multi-slit spectra of 24 of the LaEG candidates with the upgraded dual beam 
Low Resolution Imaging Spectrograph at the Keck I 10m telescope. Two separate multi-slit masks were used, 
with total exposure times of 10,800 s for each. We confirmed that 18 of the 24 galaxies were at the redshift of 
the radio galaxy. A histogram of the relative velocities and their spatial distribution is shown in Figure 3. The 
velocity dispersion is large (~ 824 kms -1 ) with a hint of kinematic substructure at ~ 1000 kms -1 blueward 
of the radio galaxy. A bimodal velocity distribution has also been found in the z = 2.156 HzRG proto-cluster 
around MRC 1138—262 (Table 2), and is fairly common even in much lower redshift galaxy clusters. The spatial 
distribution shows possible evidence for substructure as well, with small groups of 4 - 6 LaEGs to the West and 
South of the radio galaxy. More spectroscopic observations of the remaining LaEG candidates are planned to 
investigate this. 

3.2. TN J1338-1942 at z = 4.102 

TN J1338— 1942 was imaged using a narrow-band filter with a 60 A bandpass centered at the z = 4.102 redshifted 
Lya line (6202 A) with the ESO VLT using the same telescope and instrument set up as for MRC 0943—242. 
The total exposure time was 33,300 s, reaching a magnitude limit of NB^ B = 27.4 (3<r, 2.0" aperture) and the 
data were obtained during photometric conditions and 0.8"seeing. A broad-band B image (Bab = 27.02, 3<r, 
2" aperture) was used to identify LaEGS's in the same manner as for MRC 0943—242 . A total of 33 objects 



velocity distribution at z = Z.B1B (bin dm = 2DQ km/e) 




Figure 3. [Top Left] VLT narrow-band Lya image (grey scale; continuum not subtracted) of the z = 2.919 radio 
galaxy MRC 0943—242 with LaEG's identified (open squares). The blank area near the NW corner was occulted by a 
slit because of a bright star. [Top Right] Histogram of the velocity distribution of the 18 LaEG's + radio galaxy which 
were spectroscopically confirmed. [Bottom Right] Spatial distribution of the LaEG's (circles) + radio galaxy (square). 
Circle diameters are proportional to Lya flux and range from 0.8 and 17 xlO" 17 erg s _1 cm~ 2 . 



were identified as potential LaEG's (excluding the radio galaxy itself), 28 of which were retained after rejecting 
5 objects based on blue B — R colors (see Ref. 21 for details). 

Multi-slit spectra of 23 LaEG candidates were obtained with the ESO FORS2 spectrograph using two masks 
and ~ 33,300 s exposures each. The observations, and Occam's razor, 49 confirmed that 20 of the galaxies are 
at the redshift of TN J1338— 1942 (Fig. 4). The velocity dispersion of the galaxies (326 kms^ 1 ) is narrower 
than in the other, lower redshift proto-clusters (Table 2) and shows no evidence for substructure. The spatial 
distribution appears to be fairly random except that the radio galaxy itself, which is presumably the most 
massive galaxy, does not appear to be near the centroid of the proto-cluster. 

3.3. Conclusions for proto clusters 

Narrow-band imaging of the fields around HzRGs provides a powerful method for finding proto-clusters over 
a large range of redshifts. LaEG's represent only a small fraction (20 — 25%) of the total galaxy population in 
proto-clusters. 20 However, for faint galaxies (Iab > 25.5) at high redshifts emission-line objects such as these 
may be the only ones for which redshifts may be measurable, even with 8 - 10m class telescopes. 

The mass function of galaxy clusters and groups depends strongly on cosmological parameters and evolves 
with lookback time. 50 Observations of HzRG proto-clusters will be very helpful in investigating cosmological 
models in that they may be used to determine the masses of galaxy clusters and groups at redshifts which can 
not be reached by other methods (X-ray selection, gravitational lensing) , or which are prohibitive in observing 
time ('blind' multi-color surveys with no prior redshift information). 




Figure 4. [Top Left] VLT narrow-band Lya image 21 (grey scale; 6.4' x 6.2') of TN J1338-1942 with 20 spectroscopically 
confirmed LaEG candidates identified (circles). [Top Right] Histogram of the velocity distribution of the LaEG's. [Bottom 
Right] Spatial distribution of the LaEG's (circles) + radio galaxy (square). Circle diameters are proportional to Lya 
flux and range from 0.8 and 17 xlO -17 erg s _1 cm -2 . 



The over-densities in the HzRG proto-clusters are similar to that in the SSA22a over-dense region and lead 
to similar estimates of the cluster masses of 10 14 — 10 15 M Q . In fact, the luminosity functions and lifetimes of 
luminous radio sources are consistent with every such over-dense field hosting a massive galaxy that will be an 
active radio source at least once in its lifetime. 21 



4. FUTURE WORK 

Large 8 - 10m telescopes for the first time unlock the potential of narrow-band imaging as a cosmological tool. 
Giant emission-line nebulae are now known to exist out to z ~ 4, some with obvious ionizing sources, some 
without (Table 1). By studying their properties we may learn about the formation of galaxies, the feedback from 
AGN/starburst winds, and the chemical enrichment of intra-cluster media. Radio-loud galaxies seem endowed 
by rather spectacular nebulae, probably they are the most massive forming galaxies, but other types of objects 
need to be studied as well. For example, if quasars would become active during the very early stages of galaxy 
formation, then one might expect them to boost the Lya emission in the surrounding gas, which might be used 
to learn more about the galaxy formation process. 51 

Emission-line imaging of the fields around targeted high redshift objects can push searches for galaxy proto- 
clusters to much higher redshifts than would otherwise be possible. Again radio galaxies are ideally suited for 
this because they live in the most massive galaxies which live in the most over-dense regions. 

To fully exploit emission-line imaging as a cosmological tool will require large field-of-view (0.5 x 0.5 sq.deg) 
tunable filters, 52 Imaging Fourier Transform Spectrometers, 53 or other types of multi-wavelength '3-D' imaging 
devices. 54 Large (> 200 cm 2 ) custom made narrow-band filters are expensive, and it is an unfortunate fact 



that the largest telescopes also require the largest filters, even for modest sized fields of view (273 cm 2 for a 
81 sq.arcmin FOV in the case of DEIMOS at Keck for example). This problem will be even more acute for 
the next generation of 'extremely' large telescopes such as the CELT 30m. It is also clear that the selection of 
emission-line objects provides only partial insight into the galaxy populations and kinematics of proto-clusters, 
and obtaining spectroscopic absorption line redshifts of faint galaxies will be very difficult. The most powerful 
observational tool here would again be a 3-D imaging device with selectable, medium sized (few 100A wide) 
bandpasses so that accurate photometric redshift measurements can be made. 
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